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different  boundary  layer  effects,  and  these,  in  turn,  are  iden- 
tical with  the  changes  likely  to  be  produced  by  polymer  injections., 
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scope  of  the  present  study.  However,  it  is  argued  that  this 
interaction  must  necessarily  be  dependent  on  the  visco-elastic 
properties  of  the  polymer.  For  example,  it  is  shown  that  when 
the  observed  changes  in  the  lift  force  are  normalized  with 
respect  to  the  local  velocity  at  the  location  of  the  injection 
slit,  the  data  for  different  polymers,  different  injection 
velocities  and  different  angles  of  attack  can  all  be  correlated 
in  terms  of  the  single  nonditnensional  parameter  Vj  r /o  where  Vj 
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ABSTRACT 

Existing  data  on  the  effects  on  lift  of  symmetric  hydro- 
foils due  to  the  injection  of  polymer  solutions  on  to  their 
surfaces  show  that  the  lift  can  either  decrease  or  increase 
depending  on  the  polymer,  injection  velocity,  location  of  in- 
jection, side  of  injection  and  angle  of  attack.  A unified, 
albeit  semierapirical,  approach  is  suggested  in  the  present  re- 
port for  correlating  and  explaining  the  apparently  contradic- 
tory data.  It  is  proposed  that  the  observed  lift  changes  may- 
be due  to  the  influence  of  the  injections  on  the  boundary- 
layer  displacement  effect  that  is  known  to  reduce  the  circu- 
lation around  hydrofoils. 

Four  separate  boundary-layer  effects  due  to  injection  are 
identified  and  discussed,  namely  those  due  to  changes  in  the 
effective  angle  of  attack,  the  effective  camber,  the  circula- 
tion (in  addition  to  the  changes  in  circulation  associated 
with  the  previous  two  effects),  and  the  thickness  distribution 
Classical  methods  of  airfoil  theory  are  utilised  to  calculate 
the  “effective"  hydrofoil  shape  that  would  have  produced  the 
observed  pressure  distribu  Ions.  This  effective  shape  is 
utilised  to  identify  the  different  boundary  layer  effects,  and 
these,  in  turn,  are  identical -with  the  changes  likely  to  be 
produced  by  polymer  injections. 
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No  explicit  treatment  of  the  exact  nature  of  the  inter- 
action between  the  injected  polymer  flow  and  the  boundary  layer 
or  the  external  pressure  gradient  is  given,  since  this  is  out- 
side the  scope  of  the  present  study.  However,  it  is  argued  that 
this  interaction  must  necessarily  be  dependent  on  the  visco- 


elastic properties  of  the  polymer.  For  example.  It  is  shown 
that  when  the  observed  changes  in  the  lift  force  are  normalised 
with  respect  to  the  local  velocity  at  the  location  of  the  'in- 
jection slit,  the  data  for  different  polymers,  different  in- 
jection velocities  and  different  angles  of  attack  car.  all  be 
correlated  in  terms  of  the  single  nondimensional  parameter 
VjT/e  where  is  the  injection  velocity,  t is  the  polymer  re- 
laxation time  and  o is  the hydrof oil  chord.  A similar  Corre- 
lation is  also  given  for  the  changes  in  drag  due,. to  injection. 
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:i  FOREWORD 

■i  . 

The  present  report  summarizes  the  work  that  has  been  carried 
out  over  the  last  four  years  at  HYDRONAUTICS*  Incorporated*  under 
the  support  of  the  Office  of  Naval  Research*  on  the  lift  and  drag 
effects  that  arise  due  to  the  injection  of  polymer  additives  on 
to  the  surface  of  hydrofoils.  The  work  has  been  primarily  ex- 
perimental and  an  effort  is  made  herein  to  offer  a unified*  albeit 
semiempirical,  interpretation  to  the  considerable  body  of  data 
that  has  been  acquired. 

Much  of  the  early  work  on  lift  and  drag  effects  were  carried 
out  under  the  general  direction  of  Mr.  Marshall  P.  Tulin*  and  a 
number  of  the  correlations  given  in  Section  IV  are  based  on  theo- 
retical concepts  originated  by  him.  Most  of  the  pioneering  effort 
at  HYDRONAUTICS*  Incorporated  on  lift  effects  on  hydrofoils  were 
carried  out  by  Dr.  Daniel  H.  Fruman. 

Technical  monitoring  for  the  program  was  provided  by  Messrs, 
Stanley  Doroff  and  Ralph  Cooper  of  the  Office  of  Naval  Research. 
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I.  INTRODUCTION 

It  is  well  known  that  when  dilute  polymer  additive;-  are 
present  In  the  external  flows  over  two-dimensional  hydrofoils 
they  produce  changes  In  both  the  drag  and  lift  characteristic:! 
of  hydrofoils.  However,  while  the  drag  reduction  achieved  Is 
generally  accepted  as  being  due  to  a thickening  of  the  viscous 
sublayer  adjacent  to  the  wall,  the  reasons  for  the  lift  changes- 
are  far. from  clear.  Indeed  the  available  experimental  data  are 
themselves  >\ul  t-c  ambiguous  and  contradictory  as  to  the  exact 
nature  of  the  lift,  effect. 

Since  Wu's  early  finding  on  the  possibility  of  a reduction 
of  (pump)  Impeller  thrust  In  additive  solutions L , related  stud- 
ies with  propellers  and  hydrofoils  have  been  reported  by 

' > • I i, 

Kowalski  , Wolff  and  dahn  , l.ehtnan  and  Duessmann  , i’arpkaya  , 
l*’rumau,  Dundarnm  and  Daugard  , Kmman,  Tulin  and  l.iu1  , and 
Dinnarwalla  and  bunds ram'1  . f rom  d l r e o l-  measurements, 

Kowalski*  found  that-  polymer  additives  decreased  the  thrust  and 
increased  the  torque  of  a propeller  reuniting  In  a decreased 
efficiency.  On  the  other  hand,  for  a hydrofoil  with  polymer 
ejection,  Wolff  and  dahn'  found  a lift  reduction,  whereas  ; china u 
and  Suessmanu'  reported  that  lift  either  Increases  or  decreases 
depending  upon  whether  the  polymer  is  ejected  from  the  suction 
side  or  from  the  pressure  side. 

Rlnee  l>k'\,  UYPKONAUTU’D,  Incorporated  has  been  cn trying 
out-  a research  program,  under  the  sponsorship  of  the  Off  lee  of 
Naval  Keseareh,  on  various  aspects  of  the  effects  of  polymer 
additives  on  fluid  flows.  During  the  last  four  years.,  the  efforts 
have  centered  around  Investigations  of  the  changes  In  lift,  and 
drag  that,  occur  on  hydrofoils  when  polymer  additives  are  injected 
tangentially  Into  the  flow  from  various  locations  on  the  foil  sur- 
face. Up  to  dale  the  Investigations  have  been  predominantly 
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experimental,  with  the  tests  having  been  conducted  on  10-  and 
20-cm  chord,  symmetric  hydrofoils  in  the  HYDRONAUTICS  High  Speed 
Channel.  Measurements  have  included  drag  and  lift  forces,  as 
well  as  pressure  distributions  with  and  without  injections.  The 
experimental  studies  have  yielded  a large  body  of  valuable  data 
which  can  assist  in  the  understanding  of  the  effects  of  polymer 
additives  on  external  flows.  The  objective  of  the  research  de- 
scribed in  the  present  report  is  to  utilise  this  body  of  data, 
as  well  as  other  data  that  is  available  in  the  literature,  to 
attempt  to  develop  a theoretical  framework  which  can  be  used  to 
analyse  and  Interpret  the  observed  effects. 

In  order  to  set  the  objectives  of  the  present  report.  In  a 
proper  perspective,  a brief  discussion  of  the  experimental  ob- 
servations, with  emphasis  on  the  most  recent  results,  :ls  first 
given  In  Section  11,  followed  by  a detailed  description  of  ob- 
jectives. A discussion  of  the  various  physical,  mechanisms  that 
have  been  postulated  In  the  literature  to  explain  tin'  observed 
lift  effects,  arc  discussed  In  Section  111  In  the  light  of  the 
most  recent  experimental  data.  The  mechanism  that  is  considered 
to  be  the  most  plausible,  as  well  as  the  reasons  for  considering 
it  to  be  so,  are  also  described  in  Section  111,  A semi -emp l rl on  I 
analysis  is  then  given  In  Section  IV  and  Is  utilised  to  correlate 
the  data.  Finally,  some  concluding  remarks  are  given  in  Section 
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II.  FEATURES  OF  THE  EXPERIMENTAL  RESULTS  AND  OBJECTIVES 

OF  THE  STUDY 

Most  of  the  experimental  studies  conducted  at  HYDRONAUTICS, 
Incorporated  on  the  lift  effects  of  polymer  additives  were  carried 
out  with  a 10-cm  chord  NACA  63* -020,  two-dimensional  hydrofoil 
placed  in  the  HYDRONAUTICS'  High  Speed  Channel  (HSC).  In  the 
early  tests,  the  forces  on  the  model  were  measured,  by  the  use  of 
block  gauges,  for  the  cases  with  and  without  injection  of  Poiyox 
WSR301  as  well  as  water  (so  as  to  provide  a "baseline"  condition 

for  assessing  the  effects  of  the  polymer  alone).  These  early 

6 7 v 

studies  indicated  that,  under  certain  conditions,  polymer  in- 
jection not  only  leads  to  a drag  decrease  but  also  to  a lift  in- 
crease, so  that  significant  increases  in  the  lift/drag  ratio 
could  be  realised.  In  general,  it  was  found  that  a lift  increase 
accompanied  a drag  decrease  when  the  polymer  Injection  was  on  the 
suction  side  of  the  hydrofoil,  and  that  both  drag  and  lift  de- 
creased when  the  polymer  injection  was.  on  the  pressure  side  of 
the  hydrofoil. 

In  an  attempt  to  understand  the  phenomenology  responsible 
for  the  observed  lift  behavior,  tests'  were  undertaken  to  measure 
the  pressure  distributions  on  the  hydrofoil  for  the  eases  with 
and  without  polymer  injection.  These  preliminary  tests  indicated 
that  polymer  injection  on  one  surface  (suction  or  pressure)  of 
the  hydrofoil  changes  the  pressure  distribution  on  bo t.h  surfaces  , 
though  the  effect  is  more  pronounced  on  the  surface  on  which  the 
injection  is  made.  It  was  also  found  that  the  pressure  distri- 
butions were  consistent  with  the  lift  forces  measured  by  the  block 
gauges;  that  is,  the  forces  obtained  by  Integrating  the  measured 
pressure  distribution  were  the  same,  within  accepted  experimental 
error,  as  those  measured  by  the  block  gauges. 
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Based  on  these  experimental  results,  Fruman , Tulin  ar.d  Liu*7 
concluded  that  the  observed  lift  increases  could  not  bo  explained 
in  terms  of  changes  in  the  boundary-layer-separation  point,  since 
the  pressure  distributions  revealed  that  injection  did  not  signifi- 
cantly alter  the  pressure  distributions  in  the  trailing-edge  region. 
These  authors  argued  that  the.  observed  lift  effect  may  not  be  re- 
lated directly  to  the  drag-reduction  phenomena.  There  were  also 
indications  from  the  trends  in  the  data*  that  the  observed  lift 
effect  may  be  strongly  dependent  on  the  viscoelastic  behavior  of 
the  polymer.  One  direct  method  of  examining  whether  viscoelastic 
effects  are  responsible  for  the  observed  lift  changes  is  to  con- 
duct tests  with  different  polymers  of  differing  viscoelastic 
behavior  under  otherwise  identical  experimental  conditions;  the 
results  of  such  tests  are  described  in  Reference  8 for  three 
different  polymers,  namely,  Polyox,  Polyacrylamide  and  Jaguar. 

Since  these  results  formed  the  basis  for  the  analyses  given  in 
the  present  report,  they  are  described  briefly  below. 

As  mentioned  earlier,  all  of  the  tests  were  conducted  with 
a 10-cm  chord  symmetric  hydrofoil  and  were  performed  in  the  HYPRO- 
NAUTICS  High  Speed  Channel  (HSC)  modified  to  obtain  a two-dimen- 
sional flow.  The  hydrofoil  has  two  .OOlPY-em  wide  Injection 
slits,  one  located  at  10#  chord  and  the  other  (on  the  opposite 
side)  at  30#  chord.  Each  surface  of  the  foil  has  ton  pressure 
taps  between  l8  to  8o  percent  chord  length.  The  drag  and  lift 
curves  for  the  basic  hydrofoil  ’(without  any  Injection)  are  shown 
in  Figure  1.  The  slight  as  symmetry  it;  the  hydrofoil  caused  by 
the  presence  of  the  10-  and  30-pereent  chord  slits  on  opposite 
sides  of  the  foil  is  responsible  for  the  slight  a asymmetry  noted 
in  the  lift  behavior, 

♦These  aspects  will  be  discussed  in  more  detail  in  Section  111, 
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Figure  2 shows  a comparison  between  the  measured  pressure 
distribution  and  the  results  of  a calculation  using  potential 
theory^.  While  there  is  general  agreement  between  the  measured  •>  1 

and  calculated  pressure  distributions,  one  feature  that  is  note- 
worthy is  the  relatively  large  discrepancy  between  the  two  dis-  ? j 

tributions  near  the  trailing  edge  of  the  hydrofoil.  i .*■ 

Lift,  drag  and  pressure  measurements  were  performed  on  the 
hydrofoil  under  the  following  test  conditions: 

a.  Injected  Fluids:  - 

1.  Water  .j 

2.  200  ppm  Polyox  WSR  301  (Union  Carbide) 

3.  350  ppm  Polyacrylamide  (Polyscience  - Cat.  #2806) 

4.  1500  ppm  Jaguar  WPB  (Stein,  Hall  & Co.,  Control 
#23-0548). 

b.  Injection  Velocity  Ratios:  V^/V^;  0.1  and  0.3,  (V^  « 

11  meters/sec).  " 

c.  Angles  of  Attack;  a:  0°,  2.5°  (or  3.25°)  and  5°.  . 

d.  Injection  Sides:  Suction  (upper)  and  pressure  (lower). 

e.  Injection  Positions:  10$  and  30$  chord. 

In  all,  tests  were  conducted  for  eighty  different  combina- 
tions of  tost  conditions.  A summary  of  the  test  results  follows*: 

Drag  (Figures  3 and  4) 

Water  injection  results  in  a drag  increase,  or  at  most,  a 
slight  drag  reduction;  on  the  other  hand,  polymer  injection  always 
results  in  a drag  reduction. 

Lift  (Figures  5 and  6) 

Water  injoction  always  3eems  to  produce  a lift  force  in  a 
direction  opposite  to  the  injoction  side,  the  magnitude  of  this 

*A11  of  the  experimental  data  are  given  in  tabular  form  in  App radix  B. 


**  ■) 
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lift  change  being  relatively  larger  for  Vj/V^  = .3  and  10%  chord 
injection  cases.  Polymer  injection  produces  a lift  force  in 
either  the  same  direction  as  that  of  the  injection  or  opposite 

to  it,  depending  upon  the  polymer  (compare  Jaguar  and  Polyox  at 

0°  angle  of  attack  and  30%  chord  injection.  Figure  6),  the  rate 
of  injection  (compare  0.1  and  0.3  rates  of  injection  for  Polyox 
at  0°  angle  of  attack.  Figure  5),  the  location  of  injection  (com- 
pare 10  and  30%  chord  locations  for  Polyox,  pressure-side  injee- 

re, 

| tion  at  0.1  rate.  Figures  5 and  6),  and  the  angle  of  attack 

i (compare  0°  and  2.5°  with  5°  angle  of  attack  for  the  suction- 

1 side  injection  of  Polyox  at  V^/V*  = .3  rate,  Figure  5). 

Pressure  Distribution 

The  pressure  distribution  data  were  taken  for  all  of  the 
; different  test  conditions.  Typical  pressure  distribution  curves 

are  shown  in  Figures  7 and  8.  The  change  in  the  chordwise  pres- 
sure distribution  due  to  suction-side  Injections  of  Polyox  and 
Jaguar  at  a foil  angle  of  2. *3°  is  plotted  in  Figure  7.  The 
hydrofoil  has  pressure  taps  only  between  18  and  86%  of  its  chord 
f length j nevertheless,  the  general  trend  is  good  enough  to  make 

f the  following  observations. 

r 1.  Polyox  injection  at  V^/V^  = 0.1  results  in  a pressure 

f ' decrease  on  most  of  the  suction  side  and  in  a pressure  Increase 

on  the  pressure  sidej  hence,  one  would  expect  a lift  increase, 
j On  the  other  hand.  Polyox  injection  at  V^/V^  = 0.3  results  in  a 

pressure  increase  on  most  of  the  suction  side  and  in  a pressure 
decrease  on. the  pressure  sidej  hence,  one  would  expect  a lift  de- 
crease. These  observations  are  consistent  with  the  results  given 
in  Figure  5. 


2.  Jaguar  injection  at  V^/V*  = 0.1  as  well  as  at  0.3  rato 
results  in  a pressure  decrease  on  the  suction  side  and  in  a 


£tj 

1’  , 
f _ ; 

f.  I pressure  increase  on  the  pressure  side;  hence,  a lift  increase  is 

y'-.  ;» 

i-  expected  in  both  cases.  However,  the  magnitudes  of  this  pressure 

!'  ; change  on  both  sides  are  comparatively  larger  for  ViA00  = 0.3; 

i ■;  hence,  a relatively  larger  increase  in  the  lift  is  expected  for 

?V::  that  case.  Again  these  observations  are  very  consistent  with  the 

results  given  in  Figure  5. 

[,  t-  The  occurrence  of  a sharp  negative  pressure  peak  in  the 

¥ - relative  pressure  distribution  on  the  side  of  injection  (Figures 

£ 7 and  8a)  is  a characteristic  feature  of  most  of  the  cases  con- 

sidered. Water  injection  in  some  cases  gives  rise  to  a positive 
I pressure  peak;  see  Figure  8b.  The  magnitude  ana  the  chordwise 

| extent  of  this  peak  have  an  important  bearing  on  the  net  lift 

l l change  due  to  injection. 

f ' Even  from  the  brief  description  given  above,  it  can  be  seen 

i;  . that  the  lift  behavior  is  quite  complex,  with  the  lift  change  due 

r to  injection  being  positive  or  negative  for  different  combinations 

of  conditions  involving  the  injected  polymer,  the  injection  veioc- 

£ . 

ity,  the  chordwise  location  and  side  of  the  injection,  and  so  on. 
f Also,  quite  evidently  the  exact  nature  of  the  changes  in  the 

pressure  distribution  holds  the  key  to  the  observed  differences 
in  lift  changes  for  various  combinations  of  conditions,  as  illus- 
- trated  in  Figures  7 and  8.  It  is  clear  that  any  theoretical 

explanation  for  the  lift  changes  will  have  to  account  not  only  for 
?■  ; the  rather  striking  changes  in  the  pressure  distribution  due  to 

polymer  injection,  ait  also  for  the  changes  due  to  water  injection, 
l \ The  development  of  a comprehensive  theory  to  encompass  all  of  the 

factors  involved  is  undoubtedly  quite  a complex  task,  and  this  is 
; not  the  objective  of  the  present  study.  Rather,  the  objective 

here  is  to  develop  a semi-empirical,  phenomenological  approach 
useful  for  correlating  the  observed  features  of  the  data,  appeal- 
ing where  necessary  to  dimensional  reasoning. 


; 
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In  the  next  section,  the  various  mechanisms  that  have  been 
proposed  in  the  literature  as  explanations  for  the  observed  lift 
effects  are  considered,  with  emphasis  on  the  one  we  consider  most 
plausible  based  on  the  trends  in  the  data  described  herein. 
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III.  PLAUSIBLE  MECHANISMS  FOR  LIFT  CHANGES 


As  mentioned  earlier,  based  on  their  observations,  Fruman, 
Tulin  and  Liu'  examined  several  possible  explanations  for  the  ob- 
served behavior,  A change  in  boundary-layer  separation  behavior 
was  discounted  as  a feasible  explanation,  since  the  measured  pres- 
sure distributions  do  not  display  any  evidence  of  this  (the  major 
pressure  changes  occur  somewhat  aft  of  the  slit  rather  than 
near  the  trailing  edge).  These  authors  then  considered  boundary- 
layer  thinning  on  the  side  on  which  the  injection  is  made  as  a 
possible  explanation.  The  decreased  boundary-layer  thickness  at 
the  trailing  edse  of  the  hydrofoil  can  bo  directly  related  to  the 
drag  reduction,  since  the  change  in  the  momentum  flux  in  the  wake 
must  equal  the  change  in  drag.  In  turn,  the  reduced  thickness  of 
the  boundary  layer  at  the  trailing  edge  on  the  side  on  which  the 
injection  is  made  can  be  related  to  a change  In  the  effective 
angle  of  attack.  It  is  relevant  to  note  that,  according  to  the 
above  explanation.  Injection  on  the  pressure  side  will  lead  to  a 
decrease  in  the  trailing-edge  boundary-layer  thickness  on  this 
side  and  a consequent  dec rease  in  the  effective  angle  of  attack, 
whereas  injection  on  the  suction  side  will  lead  to  the  opposite 
effect— a result  which  is  consistent  with  observed  lift  behavior 
in  most  (but  not  all)  of  the  cases. 

Fruman,  Tulin  and  Liu  calculated  the  reduction  in  trailing- 
edge  boundary- lay or  thickness  using  the  measured  drag  reduction 
and  calculated  the  change  in  lift  corresponding  to  the  change  in 
the  effective  angle  of  attack  by  using  the  measured  lift-curve 
slope.  They  found  that  in  all  but  a few  cases,  the  measured  lift 
changes  were  considerably  larger  than  those  that  can  be  predicted 
using  the  mechanism  described  above.  As  will  be  seen  later,  a 
change  in  the  effective  angle  of  attack  due  to  boundary- layer  thin- 
ning at  the  trailing  edge  is  only  a part  of  the  total  effect  on  the 
boundary  layer  of  the  polymer  injection. 
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Fruman,  Tulin  and  Liu*^  then  considered  viscoelastic  effects 
as  possible  explanations  for  the  observed  changes.  They  pointed 
out  that  when  the  observed  lift  forces  are  plotted  against  the 
logarithm  of  the  free  stream  velocity,  a straight-line  behavior 
results  in  a manner  analogous  to  the  behavior  noted  in  pitot-tube 
measurements  in  flows  containing  polymer  additives.  Such  a plot 
also  indicates  that  there  is  a threshold  velocity  below  which  the 
lift  effect  does  not  appear,  the  actual  value  of  this  threshold 
velocity  being  dependent  on  the  test  conditions.  However,  when 
the  results  are  plotted  in  terms  of  the  local  velocity  at  the 
injection  slit,  a single  threshold  velocity  results.  These  ob- 
servations lend  credence  to  the  concept  that  viscoelastic  effects 
may  be  responsible  for  the  observed  lift  effects.  Moreover,  the 
fact  that  under  otherwise  identical  experimental  conditions  differ- 
ent polymers  lead  to  significantly  different  lift  changes,  also 
lends  credence  to  the  concept  that  viscoelastic  effects  are 
responsible,  at  least  in  part,  for  the  observed  effects, 

l ^ 

The  anamolous  lift  effect  was  also  noted  by  farpkaya"  In  his 
tests  on  two-dimensional  hydrofoils  immersed  In  dilute  homogeneous 
solutions  of  polyethelene  oxide,  and  he  concluded  that  the  lift- 
to-drag  ratio  of  the  hydrofoils  may  either  increase  or  decrease 
depending  on  their  stall  characteristics,  garpkaya  fhoorl.ed 
that  the  observed  effect  may  be  due  to  the  different  influences 
of  the  polymer  on  the  boundary -layer  character Is  ties  on  the  top 
and  bottom  surfaces  of  the  hydrofoil  and  the  consequent  change  in 
the  circulation  defect  (that  is,  the  difference  between  the  actual 
circulation  and  that  predicted  by  the  Kutta  condition)  known  to 
exist  in  real  flows. 

The  hypothesis  that  the  observed  effects  may  be  due  to 
boundary-layer  interaction  becomes  quite  plausible  when  one  com- 
pares the  actual  measured  pressure  distributions  under  undisturbed 
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conditions  (that  is,  in  the  absence  of  injection)  with  theoreti- 
cal computations  based  on  thin-airfoil  theory;  see  Figure  2.  It 
can  be  seen  that  there  are  significant  differences  between  the 
calculated  and  measured  distributions,  especially  near  the  trail- 
ing edge  of  the  hydrofoil;  these  differences  are  typical  of  such 
comparisons,  are  well  known  in  the  literature  (see  References  10- 
13,  for  example),  and  are  attributed  to  boundary-layer  effects. 

The  potential-flow  streamlines  around  the  hydrofoil  are  displaced 
outward  not  only  due  to  the  thickness  distribution  of  the  hydro- 
foil, but  also  due  to  the  boundary  layer  on  the  foil  surface. 

Hence,  better  results  would  be  obtained  if,  in  the  computations 
instead  of  using  the  actual  hydrofoil  shape,  an  altered  shape 
in  which  the  boundary-layer  displacement  thickness  is  added  to  the 
shape  is  used. 

In  general,  the  boundary  layer  on  the  upper  surface  of  the 
hydrofoil  will  grow  faster  than  that  on  the  lower  surface  because 
of  the  adverse  pressure  gradient  existing  on  the  top  surface  from 
the  minimum-pressure  point  onward.  Therefore,  the  effective  hydro- 
foil shape  (that  is,  the  actual  shape  plus  the  displacement  thick- 
ness) will  have  a slightly  turned  up  trailing  edge  at  a relatively 
small  positive  angle  of  attack.  The  differences  between  the  com- 
putational results  and  the  measurement  near  the  trailing  edge  of 
the  foil  are  directly  attributable  to  the  effect  mentioned  above 
and,  indeed,  excellent  agreement  between  the  two  results  are  ob- 
tained if  boundary -lay or  displacement  effects  are  included  In  the 
oompu ta  t ion 1 1 9 1 * . 

The  relevance  of  the  above  remarks  in  the  present  context  is 
that,  as  can  be  seen  from  Figure  2,  the  boundary- layer  displacement 
effects  can  have  a significant  effect  on  the  pressure  distributions, 
these  effects  being  considerably  larger  than  the  observed  lift  ef- 
fects caused  by  polymer  injections.  Thus,  even  small  changes  in  the 
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dynamic  evolution  of  the  boundary  layer  on  the  surface  on  which 
the  injection  is  made  can  be  expected  to  produce  changes  in  the 
pressure  distributions  of  the  type  observed  in  the  experiments. 

One  other  important  feature  of  the  relative  pressure  dis- 
tributions which  also  indicates* that  the  observed  effect  may  be 
due  to  boundary-layer  interaction  is  the  relatively  sharp  de- 
crease in  the  pressure  coefficient  some  distance  aft  of  the  injec- 
tion position.  Tills  relatively  sharp  negative  peak  in  the  pres- 
sure distribution  is  a characteristic  feature  of  most  of  the 
cases  considered,  and  always  occurs  o*  the  same  side  on  which  ;he 
injection,  is  made  (see  Figures  7 and  8).  As  already  noted,  and 
as  can  be  seen  from  Figures  7 and  3,  the  magnitude  and  ehordwise 
extent  of  this  peak  have  an  important  bearing  on  the  magnitude  as 
well  as  the  direction  of  the  net  lift  change  that  results  due  to 
the  polymer  injection. 

In  many  eases,  as  in  the  case  of  Poly ox  injection  illustrated 
in  Figure  7,  there  are  positive  pressure  regions  before  and  fe l low- 
ing the  negative  pressure  peaks.  Also,  the  pressure  distribution 
is  affected  everywhere  on  the  foil  surface  regardless  of  the.  side 
or  location  at  which  infection  is  made.  Moreover,  the  actual 
location  at  which  infection  is  made  seems  to  have  little  erne 
influence  on  the  location  of  the  pressure  peak,  the  latter  appar- 
ently being  ImMueneed  more  by  the  has  I e undisturbed  rrossnre  dis- 
tribution on  the  hydrofoil. 

It  is  important  to  note  that  water  injection  also  leads  to 
peaks  in  the  pressure  distribution,  us  illustrated  In  Figure  3b. 
litis  figure  shows  the  changes  in  the  pressure  distributions  re- 
sulting from  water  injection  at  an  angle  of  attack  of  0°  for  lower 
side  injection  at  30y-  chord,-  For  the  corresponding  ease  of  Polyox 
injection,  there  is  a negative  peak  on  the  bottom-side  pressure 
distribution  (the  side  on  which  the  injection  is  made).  However, 

this  's  tlTHment  will  be  discussed  in  more  quantitative  terns  In 
Section  IV. 
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for  the  case  of  water  injection,  while  there  is  a small  negative 
peak  for  the  smaller  rate  of  injection,  the  pressure  peak  becomes 
positive  for  the  larger  rate  of  injection. 

The  observed  features  of  the  pressure  distributions  described 
above  are  compatible  with  the  hypothesis  that  the  lift  changes  are 
produced  by  boundary-layer  displacement  effects.  The  interaction 
between  the  injected  flow  and  the  boundary  layer  can  be  expected 
to  depend  as  much  on  the  undisturbed  pressure  distribution  on  the 
hydrofoil  as  it  does  on  the  location  of  the  injection.  For  ex- 
ample, clearly  portions  of  the  boundary  layer  in  the  favorable  and 
adverse  pressure-gradient  regions  in  the  hydrofoil  may  respond 
differently  to  the  perturbations  caused  by  the  injections. 

v 

Fruman,  Tulin  and  Liu'  have  suggested  that  the  lift  effect 
of  polymer  injection  may  be  due  to  the  fact  that  the  polymer  stream 
enters  the  flow  around  the  hydrofoil  in  the  form  of  a "swollen  jet" 
due  to  the  viscoelastic  behavior  of  polymer  solution.  While  it  is 
plausible  to  suppose  that  the  observed  peaks  in  the  pressure  dis- 
tribution may  be  due  to  these  "swollen  jets",  it  seems  unlikely  that 
the  observed  effect  is  a purely  viscoelastic  one,  since  vmior  in- 
jection also  leads  to  peaks  in  the  pressure  distributions. 

In  any  ease,  the  important  point  to  note  here  is  that  the 
local  changes  in  the  displacement  thickness  of  the  boundary  layer" 
due  to  its  interaction  with  the  injected  flow  will  lead  to  local 
changes  in  the  effective  camber  of  the  hydrofoil.  In  turn,  the 
effect  of  the  latter  changes  will  be  to  cause  localised  changes 
in  the  pressure  distribution  near  the  locations  of  maximum  inter- 
action as  well  as  to  cause  more  general  changes  -in-  the  pressure 
distribution  over  the  entire  hydrofoil  surface  due  to  changes  in 
the  circulation.  This  is  exactly  what  is  seen  in  the  data. 


«JV>^  ■ . '• 
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tY  ' As  mentioned  earlier,  based  on  his  experiments  on  hydrofoils 
immersed  in  dilute,  homogeneous  polymer  solutions,  Sarpkaya^  has 
■suggested  that  the  observed  lift  increases  may  be  due  to  the 
restoration,  by  polymer-induced  changes,  of  some  of  the  circula- 
tion reduction  created  by  boundary-layer  displacement  effects. 

The  present  results  on  the  effects  of  the  polymer  v/hen  it  is  in- 
jected at  the  hydrofoil  surface  also  seems  to  suggest  that  the 
■observed  lift  changes  may  be  due  to  a boundary -layer  displacement 
‘effect  caused  by  the  injection,  especially  in  view  of  the  differ- 
ences between  the  actual  measured  pressure  distributions  under 
^undisturbed  conditions  (that  is,  in  the  absence  of  injection)  and 
theoretical  computations  based  on  potential-flow  theory 9;  sec  Figure 
2.  However,  the  observed  effect  is  also  evidently  a function  of 
the  characteristic  relaxation  time,  r , of  the  polymer,  or,  more 
appropriately,  the  non-dimensional  parameter  VjT/t,  where  Vj  is  the 
injection  velocity  and  i is  a characteristic  hydrofoil  dimension. 
The  data  also  show  that  the  boundary-layer  displacement  effect  ..must 
also  be  a function  of  the  side,  location  amt  velocity  of  the  polymer 
injection. 

The  hypothesis  given  above  can  be  verified  directly,  since  the'., 
measured  pressure  distributions  can  be  analysed  using  classical 
airfoil  theory,  and  the  "effective"  hydrofoil  shape  that  will 
produce  the  observed  pressure  distributions  can  be  calculated. 

This  effective  hydrofoil  shape  can  then  be  viewed  in  terms  of  a:. 
change  in  the  evolution  of  the  boundary  layer  around  the  actual 
hydrofoil.  By  these  means,  correlations  can  be  sought  between  the 
observed  effects  and  the  test  variables  in  terms  of  information 
that  exist  in  the  literature  on  the  behavior  of  boundary  layers 
under  favorable  and  adverse  pressure  gradients,  and  under  the  in- 
fluence of  injections  and  various  surface  perturbations.  This  is 
the  approach  that  is  utilised  in  the  next  section. 
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IV.  SEMI -EMPIRICAL  CORRELATIONS  FOR  THE  OBSERVED  EFFECTS 
IV. 1 Pressure  Distribution  and  Lift  Effects 


Kw 


It  was  suggested  in  the  previous  section  that  the  observed 
lift  effects  are  due  to  a change, caused  by  the  injection  of  water 
or  a polymer  at  the  surface  of  the  hydrofoil,  in  the  boundary- 
layer-displacement  phenomenon  which  causes  the  observed  lift  in 
real  flows  to  be  less  than  that  predicted  by  potential  theory. 

In  essence,  the  faster  growth  of  the.  boundary  layer  on  the  suc- 
tion side  of  the  hydrofoil  compared  to  that  on  the  pressure  side 
leads  to  an  effective  hydrofoil  shape  which  lias  a slightly  turned-up 
trailing  edge.  In  turn,  this  leads,  to  a reduction  in  effective 
camber  and  hence  also  in  circulation  and  in  lift. 

I A Q 

It  is  well  Known  97  that  reasonable; -agreement  can  be  ob- 
tained between  measured  pressure  distributions  on  a hydrofoil 
and  the,  predictions  of  potential  theory,  except  near  the  trail- 
zing  edge,  when  the  measured  lvalues  of  the  lift  coefficient, 
rather  than  the  theoretical  values,  are  used  to  normalise  the 
pressure  coefficients. - This  is  the  basis  for  the  comparison 
shown  in  Figure  f!.  Fomowhat  better  agreement  between  the  mea- 
sured ami  calculated  pressure  distributions  can  be  obtained  by 
recognising  that  the  sinaller-tliau-potential  lift  that  is  mea- 
sured must  necessarily  Imply  that  the  actual  circulation  around 
the  hydrofoil  is  smaller  than  that  predicted  by  the  Kutta  con- 
dition, Hence,  the  distribution  of  vortices  used  to  represent 
the  hydrofoil  must  correspond  to  the  reduced  circulation  rather 
than  to  that  predicted  by  the  Kutta  condition. 


r - rk  (i  - k) 


Ci3* 


liquation  tlJ  is  equivalent  to  Cj4  « Cm  (1  - k)  where  CI(  and 
Cp  are  respectively,  the  measured  ana  predicted  values  of 

thfe  lift  coefficient. 
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h = 1 - I k 


where  F and  F^  are,  respectively,  the  actual  value  of  the 
circulation  and  that  predicted  by  the  Kutta  condition. 

It  can  be  shown-^36,12 that,  for  a symmetric  hydrofoil, 
the  induced  velocity  distribution  due  to  the  modified  vor- 
tieity  distribution  is  given  by. 


l+(da/dx)y 


'x  1-2 


1 + I3(x)  -kl4  (x) 


m 

} • 


-C33* 


where  V is  the  free-stream  velocity,  a is  the  angle  of 


attack  and  the  integrals.  Ia  and  I*  are  given  by, 

i 


mm  -i  J[S  - 


14  (ko  ) « t4 


H 


z dx 
u-x  r> 


In  the  above  equations  z - s(x)  represents  the  thickness 
distribution  of  the  hydrofoil  and  all  the  lengths  have  been 
normalised  through  division  by  the  hydrofoil  chord,  c.  in 
Equation  l$]  the  positive  sign  applies  to  the  upper  surface 
of  the  hydrofoil  while  the  negative  sign  applies  to  the  lower 
side. 


In  Equation  C'3J , the  parameter  y,Li(ds/dx)a  represents 
the  ratio  between  the  (predicted)  velocity  on  the  chord line 
and  the  actual  velocity  an  the  foil  surface;  its  use  avoids, 
for  a round-nosed  hydrofoil,  the  occurrence  of  an  infinite 
s ue  tion  velocity  at  the  leading  edge^’K  Since  Equation  [3] 


. p • ' • Mo 


m 


efhTs  equation  is  found  by  representing  the  hydrofoil  surface 
by  appropriate  distribution  of  vortices  and  sources  (see  Appendix  A). 
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violates  the  Kutta  condition*  it  predicts  an  Infinite  velocity 
at  the  trailing  edge;  however*  this  can  be  eliminated  if*  as 
suggested  by  Preston10  and  Spence1^*  the  velocities  at  the  edges 
of  the  boundary  layer  along  normals  to  the  upper  and  lower  sur- 
faces at  the  trailing  edge  are  used*  along  with  a judicious  fair- 
ing procedure*  instead  of  the  velocity  at  -the  trailing  edge 
' itself. 

Equation  [3]  represents  the  velocity  distribution  due  to 
the  effect  of  angle  of  attack*  and*  for  a symmetric  hydrofoil* 
the  velocity  due  to  thickness  distribution  has  to  be  added  to 
this  to  obtain  the  total  velocity  distribution.  The. latter  veloc- 
ity distribution  can  be  written  as 


V = 


Vco  (l+lj  ) 

rl+(dz/dx)“' 


where 


(*>  ) 


(ds/dx)  dx  ^ 
xj-x  dx 


The  velocity  due  to  thickness  distribution  is  also 
tabulated  in  Reference  9 for  a number  of  NACA  foil  shapes 
in  general*  and  for  the  specific  hydrofoil  under  consideration 
here  in  particular. 

Equations  [3]  and  [6]  can  be  used*  along  with  the  simple 
relation, 

-‘-ter- 

to  calculate  the  pressure  coefficient*  C * on  either  side 
of  a symmetric  hydrofoil  at  a specified  angle  of  attack  and 
at  a arbitrary  value  of  the  circulation  (that  is*  at  a value 


If 
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of  the  circulation  that  does  not  necessarily  correspond  to  that 

given  by  the  idealized  Kutta-Joukowski  condition).  The  oaleu- , 
lated  pressure  distributions  for  the  specific  hydrofoil  shape 
utilized  in  the  present  study  are  shown  in  Figures  9 and  10  for 
an  angle  of  attack  of  5°  and  for  a value  of  the  circulation 
which  corresponds  to  the  measured  value  of  the  lift  coefficient. 

In  the  approximate  calculation,  the  integrals  given  in.  Equations 
[4]  and  [5]  were  not  calculated j rather,  the  expression  within 
the  flower  brackets  in  Equation  [ 31  was  equated  to  unity  since, 
this  is  sufficiently  accurate  for  the  present  purpose*.  Also, 
the  tabular  values  given  in  Reference  9,  rather  than  Equation 
[6],  were  used  to  calculate  the  induced  velocity  due  to  the 
thickness  effect. 

A comparison  between  the  measured  and  calculated  values  of 
the  differences  in  the  pressure  coefficients  on  the  upper  and 
lower  surfaces  of  the  hydrofoil  is  shown  in  Figure  9.<  while  the  ■ 
pressure  coefficients  themselves  are  shown  in  Figure  10.  Values 
corresponding  to  the  Kutta-Joukowski  value  of  the  circulation 
are  also  shown  in  these  figures  for  purposes  of  comparison.  It 
ean  be  seen  that  the  calculations  for  the  reduced-circulation 
case  give  very  good  agreement  with  the  measurements,  Including 
in  the  region  of  the  trailing  edge,  which  is  noteworthy  in  view 
of  a quite  approximate  method  of  calculation. 

The  main  point  to  be  derived  from  the  comparisons  shown  in 
Figures  9 and  10  Is  that  while  the  considerably  smaller  lift  that 
is  measured  in  the  tests  (in  the  absence  of  injection)  may  be  caused, 
in  part,  by  a lack  of  two-dimensionality  in  the  flow  or  by  wall  inter- 
ference effects,  boundary -layer  displacement  effects  may  also  be  of 
importance,  hence,  it  is  plausible  that  even  if  polymer  Injection 
caus.es  only  small  changes  In  the  boundary -layer  thickness,  the 
observed  magnitudes  of  lift  and  pressure  changes  can  result. 

*Th:Ls  approximation  is  equivalent  to  making  the  thi n-airfo l 1 assump- 
tion when  determining  the  vortex  distribution  on  the  chord lino j see 
Appendix  A. 
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Estimates  of  the  types  of  changes  that  will  occur  in  the 
.pressure  distributions  on  a symmetric  hydrofoil  due  to  the  ob- 
served lift  changes  caused  by  injection  can  be  made  using  Equa- 
• tion  [33.  Again,  as  before,  making  the  approximation  that  the 
terms  within  the  flower  brackets  in  Equation  [3]  are  equal  to 
unity,  this  equation  and  Equations  [l],  [2]  and  [8]  yield  an 
expression  for  the  change,  AC^  , in  the  pressure  coefficient 
- due  to  a change,  AC^  , ,in  the  lif-t  coefficient,  as 

(VA»)  acl 

AC  = q.  —=  =======■ -a=s=======  • — j [93 

p /l+(dz,/dx)2  • /x(l-x)  LTT 


where,  now,  the  negative  sign  applies  to  the  top  surface  and 
the  positive  sign  to  the  lower  surface. 

In  Equation  [93 j V represents  the  velocity  on  the  upper 
or  lower  surfaces  of  the  hydrofoil  in  the  absence  of  injection 
and  can  be  computed  using  the  measured  pressure  distribution 
in  the  absence  of  inj.ecti.on;  these  basic  pressure  distributions 
for  the  four  angles  of  attack  considered  in  the  present  study  . 
are  shown  in  Figure  11. 

Comparisons  between  the  measured  changes  in  pressure  dis- 
tributions duo  to  injection  and  those  calculated  using  Equa- 
tion [93  are  shown  in  Figure  12  for  two  cases.  The  first  case 
shown  is  that  for  Jaguar  injection  on  the  suction  side  of  the 
hydrofoil  at  ton  percent  chord  and  at  thirty  percent  of  1‘roo- 
stream  velocity,  with  the  angle  of  attack  being  2.9°;  as  can 
be  seen  from  Figure  5,  this  case  leads  to  a substantial  lift 
increase.  The  second  case  shown  is  that  for  water  injection, 

this  case  leads  to  a lift  decrease. 
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It  can  be  seen  from  Figure  IF  that  there  is  a very  good 
general  agreement  between  the  calculated  and  measured  distri- 
butions, except  for  the  presence  in  the,  latter  of  the  pressure 
peaks  already  discussed.  It  should  be  noted  that  in  each  of  the 
comparisons  shown  three  separate  sets  of  measurements  are  in- 
volved, namely,  the  changes  in  the  lift  force,  the  basic  pressure 
distribution  and  the  change  in  the  pressure  distribution  due  to 
injection.  Thus,  the  good  agreement  displayed  in  Figure  12  is 
evidence  of  the  self-consistency  of  the  experimental  data.  The 
fact  that  the  pressure  distributions  predicted  by  potential 
theory  using  a change  in  circulation  alone  go  a long  way  in  pre- 
- dieting  the  observed  changes  in  pressure  distributions,  appears 
to  -support  the  present  method  of  approach  of  viewing  the  lift 
and  pressure  changes  in  terms  of  changes  in  the  "effective"  shape 
of  the-  hydrofoil . 

On  the  other  hand,  when  the  observed  lift  changes  are  viewed 
merely  as  due  to  a change  in  the  effective  angle  of  attack  of  the 
hydrofoil',  the  predicted  pressure  changes  are  not  in  good  agree- 
ment with  those  observed.  This  can  be  seen  by  letting  h and  k In 
Equation  [3]  equal  to  unity  and  aero,  respectively,  so  as  to  ob- 
tain the  velocity  distribution  corresponding  to  the  ideal  value 
■ of  the  condition  around  the  hydrofoil.  Then,  the  change  in  the 
velocity  distribution  due  to  a change,  say  in  the  angle  of  attack, 
.can  be  written  as. 


AV  ® ± 


da/  dx 


[10] 


where  the  same  approximation  as  that  involved  in  Equation  [o]  hat 
been  made. 

Using  Equations  [8]  and  [lo],  the  change  in  pressure  co- 
efficients can  then  be  written  as. 
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AC  = ± 
P 


(VAJ 

fl  + (dz/dx)2 


(SCl/da)i 


[11] 


where  (SC^/Aa)m  is  the  measured  lift  curve  slope  in  the  - 
absence  of  injection.  Equation  [ll]  is  also  shown  plotted 
in  Figure  12  (in  dotted  lines);  it  can  be  seen  that,  even  the 
qualitative  behavior  of  the  measured  and  predicted  distri- 
butions are  different.  . It  should  be  noted  that  both  Equations 
[9]  and  [11]  correspond  to  the  same  change  in  circulation 
around  the  hydrofoil,  which,  in  turn,  is  deduced  from  the 
measured  value  of  the  lift  change.  However,  these  two  equa- 
tions account  for  the  changes  in  circulation  in  terms  of 
different  distribution  of  vortices  around  the  hydrofoil  chord; 
the  physical  significance  of  the  differences  between  the  two 
cases  will  be- discussed  presently. 

As  already  noted,  the  single  expression  given  in  Equation 
[9]  does  remarkably  well  in  predicting  many  of  the  observed 
features  of- the  pressure  distributions.  However , since  it  only 
accounts  for  the  changes  .in  the  pressure  distributions  due  to 
a change  in  circulation  for  a symmetric  hydrofoil,  the  equation 
fails,  not  surprisingly,  to  predict  the  observed  peak  in  the 
pressure  distribution  on  the  side  on  which  the  injection  Is 
made.  This  pressure  peak  is  evidently  associated  with  an  un- 
symmetrle  change  in  the  effective  shape  of  the  hydrofoil  and 
can  be  predicted  only  by  explicitly  accounting  for  camber 
effects.  Note  that  boundary-layer-displaeemont  effects  were 
not  accounted  for  explicitly  in  Equation  [9] > though  the  changes 
in  circulation  were  assumed  to  be  induced  by  changes  In  boundary  - 
layer  thickness. 

As  mentioned  earlier,  several  authors  have  included  the 
effects  of  boundary-layer  displacement  on  lift  explicitly.  Speci- 
fically, Preston'^0* ^ has  resolved  the  displacement  effect 
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into  symmetric  and  antisymmetric  parts  and  has  treated  each 
separately.  Preston.  '3  "change  ol'  shape"  method  is  of  consider- 
able value  in  visualising  the  effects  of  surface  injections. 
However,  before  considering  this  method  in  detail  it  is  relevant 
to  briefly  discuss  an  alternate  technique  developed  by  cpen.ee*", 
since  it  provides  some  insight  into  the  manner  in  which  boundary- 
layer  development  influences  the  outside  potential  flow.  In 
this  method  the  boundary-layer  displacement  effect  is  represented 
at  the  outset  by  an  nnsynmetric  source  distribution  and  the  cir- 
culation is  determined  by  utilising  the  condition  that  equal 
amounts  of  vortioit-y,  though  of  opposite  signs,  must  pass  into 
the  wake  from  the  boundary  layer  on  the  upper  ar.d  lower  surfaces 
of  the  hydrofoil. 

The  calcu.latier.al  store  in  fcenee 1 s method  involves  the 


transformation,  by  the  use  of  conformal  mapping,  of.  the  hydro- 
foil into  a circle;  as  such,  the  final  results  in  this  method 
involve  complex  coordinates , In  particular,  two  of  the  hey 
parameters  that  appear  in.  the  final  results  (Fount ion  -!.f  of 
Reference  13)  are  \x  and  \8  which  are,  respectively,  the  dis- 
tances in'  the  complex  plane  between  the  trailing  edge  ana  the 
edges  of  the  boundary  layers  on  the  upper  and  lower  surfaces. 
While  the  result  itself  is  too  Involved  to  consider  hers,  it  is 
relevant,  to  cite  gpenee's  observations  on  the  results. 

"This  equation  provides  a theoretical  explanation  of  cer- 
tain qualitative  boundary -layer  effects  on  lift  which  are  w\.  1 1 
recognised  experimentally.  For  an  airfoil  of  finite  trailing, 
edge  angle,  t , the  effect  of  thickening  the  upper  boundary 
layer,  keeping  the  lower  constant,  would  be  to  increase  . The 

effect  would  thus  be  to  decrease  the  circulation  and  lift, 

as  expected.  The  converse  effect  should  occur  if  X*  were  de- 
creased ( e . g . , by  Increasing  the  amount  of  laminar  flow  on  the 
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upper  surface*).  The  effect  of  increasing  and  decreasing  \z  , 
which  represents  the  lower  boundary  layer  thickness  are  not  clear 
from  the  equation,  and  would,  appear  to  depend  on  actual  rather 
than  relative  thicknesses.  It  is  expected  that  the  larger  the 
trailing-edge  angle,  the  more  marked  will  he  such  changes  in  lift 
obtainable  via  trailing-edge  thickness." 

While  the  exact  nature  of  the  interaction  between  the  in- 
jected polymer  flow  and  the  boundary  layer  is  not  known,  it  is 
clear  that  many  of  the  observations  of  Spence  quoted  above  may  be 
relevant  in  the  context  of  the  present  study.  For  example,  if 
the  effect  of  polymer  injection  on  the  upper  surface  is  imagined 
as  producing  a more  laminar  (that  is,  thinner)  boundary-layer 
behavior  on  that  surface,  then  the  result  of  the  injection  would 
be  to  yield  a lift  Increase.  This  is  what  is  observed  in  many 
cases  of  polymer  injection  studied  herein,  especially  when  the 
effects  of  polymer  injection  are  viewed  relative  to  that  of  water 
injection.  However,  it  should  be  noted  that  the  injected  flow 
will  also  interact  with  the  external  pressure  gradient,  and  that 
the  effect  of  the  strong,  adverse  pressure  gradient  present  down- 
stream of  the  point  of  injection  (see  Fig.  11)  may  be  to  increase 
the  boundary  layer  thickness.  The  net  effect  will  them  be  gov- 
erned by  the  balance  between  these  two  counteracting  influences. 

Water  injection  on  the  upper  surface  presumably  loads  to  a 
thickening  of  the  boundary  layer,  since  a reduction  in  lift  always 
accompanies  such  injection.  Also,  since  water  injection  on  the 
upper  surface  causes  a drag  increase,  it  is  difficult  to  attribute 
the  observed  results  to  any  loading-edge  effects  such  as  those 
associated  with  boundary- layer  control  through  the  use  of  loading- 
edge  slats  or  blowing^, 


^ Fmphaa i a su pp 1 i od ; not  in  the  original  citation. 
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Spence* s results  also  show  that  the  effects  on  lift  of 
changes  in  the  upper  and  lower  boundary  layer  thicknesses  are 
considerably  different.  Specifically,  changes  in  the  lower  sur- 
face boundary  layer  can  produce  either  an  increase  or  decrease 
in  lift  depending  on  the  actual  values  of  the  boundary  layer 
thicknesses  on  the  upper  and  lower  surfaces  (that  is,  on  section, 
shape,  angle  of  attack,  Reynolds  number,  etc.)  and  on  the  trail-''’" 
ing-edge  angle.  The  experimental  results  are  qualitatively  con- 
sistent with  this  observation,  since  the  changes  in  lift  due  to 
lower  surface  injection  display  considerably  fewer  trends  as 
compared  to  effects  of  upper  side  injection. 

As  mentioned  earlier,  Spence's  results  are  not  in  a form 

from  which  quantitative  estimates  can.  be  obtained  readily;  for 

1 0 1 1 

this  purpose  we  appeal  to  the  technique  developed  by  Trenton' 

Consider  a symmetric,  two-dimensional  hydrofoil  at  an  angle 
of  attack  a to  the  free  stream  as  shown  in  Figure  d^al.  hue  to 
the  more  adverse  pressure  gradients  on  the  upper  surface  of  the 
hydrofoil,  the  boundary- layer  displacement,  thickness  on  that 
surface  will  increase  at  a more  rapid  rate  than  that  on.  i he  lower 
surface.  The  effect-  of  boundary-layer  displacement-  on’  the  ex- 
ternal potential  flow  can  oe  calculated  by  adding  the  displacement, 
thickness  to  the  effective  eross-seotiona l Shane  of  the  hydro- 
foil, If  & and  5(  represent,  the  displacement-  thicknesses  of 
the  upper  and  lower  surfaces,  the  d:t splaoompnt-  effect  can  be  re- 
solved into  two  parts,  one  symmetric  and  the  other  antisymmetric , 
namely, 

* * (*U  * 6e  - * <*u  - O ' CIS] 

where  5*  represents  the  change  in  the  effective  camber  of  the 
l-  v 

hydrofoil  and  5,  the  change  in  its  thickness  distribution. 

* in  this  and  the  following  figures  the  boundary- layer-displace- 
ment effects  are  shown  highly  exaggerated  for  the  sake  of 
clarity. 
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It  is  appropriate  to  first  consider  the  effective  ea::.':.-sr 
distribution  alor.e,  since  its-  influence  on  the  lift  is  much  me r 
profound  than  that  of  the  thickness  distribution.  The  of  fee  tiv-. 
centerline  of  the  hydrofoil  is  obtained  by  adding  fcc'v  zo  its 
geometric  centerline,  and  the  effective  cross-sectional  shape  os 
the  hydrofoil  can  be  represented  as  in  Figure  13(b).  As  can.  be 
seen  from  the  figure,  the  displacement  effect  will  produce  not 
only  a reduction  in  the  effective  angle  of  attack  by  an  amount- 
equal  to, 


a = 

i 


trailing  edge 


but  also  a negative  camber  about  the  effective  chord,  AT'. 

Both  of  these  effects  will  tend  to  reduce  the  lift  produced 
by  the  hydrofoil . 

If  the  value  of  the  circulation  about  the  effective  hydro- 
foil shape  wore  governed  by  the  Kutta-Joukowsk L hypothec's,  then 
the  above  two  effects  will  be  the  only  ones  present,  however, 
due  to  the  pres  one  e of  the  boundary  layer,  the  oirou.'Uitlon  around 
the  effective  hydrofoil  is  not  given  by  tho  Kutfa-doukowsUi 
hypothesis’* ‘* , but,  rather,  has  to  bo  determined  using  ;he  con- 
dition that  equal  and  opposite  amounts  of  vortleivy  must  pass 
into  the  wake  from  the  upper  and  lower  boundary  layers,  Prcsies 
has  ishown  that  the  above  condition  is  satisfied,  to  a first  ap- 
proximation, when  tho  velocities  at  the  edges  of  the  boundary 
layer,  on  the  upper  and  lower  surfaces*  at  the  trailing,  edge  of 
the  hydrofoil  are  equal.  These  velocities  themselves  can  be 
estimated  by  using  potential  theory  with  an  arbitrary  value  for 
the  circulation  in  the  absence  of  displacement  effects,  and  the 
circulation  adjusted  to  make  the  velocities,  equal,  if  necessary, 
the  results  can  be  improved  by  iteration. 
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The  effect  of  the  symmetric  part,  6 , on  the  pressure 

distribution  can  be  calculated  by  adding  5 to  the  thickness 

□ 

distribution  of  the  hydrofoil  and  by  using  Equations  [6]  anti  [?]. 

The  effect  of  boundary-layer  displacement  on  the  lift  and 
pressure  distributions  on  a symmetric,  two-dimensional  hydrofoil 
is  thus  seen  to  consist  of  four  separate  parts.  The  first  is 
associated  with  a reduction,  in  the  effective  angle  of  attack, 
while  the  second  is  associated  with  the  n.ega tlve  camber  of  the 
effective  hydrofoil  shape.  Hie  third  effect  represents  the  re- 
duction in  circulation  as  compared  to  that  given  by  the  Kutta- 
Joukowski  criterion,  and  the  fourth  effect  represents  the  changed 
thickness  distribution.  In.  the  comparisons  shown  in  'figures  e. 
and  10,  only  the  third  of  the  above  effects  was  included  (using 
the  measured  lift)  and  the  others  were  neglected. 

The  influence  of  polymer  or  water  injection  on  the  lift  and 
pressure  distributions  can  now  be  viewed  in  terms  of  the  boundary - 
layer-displacement  effects  described  above.  Consider  the  injec- 
tion of  polymer  or  water  from  a location,  .1  , on  the  surface  of 
a symmetric,  two-dimensional  hydrofoil.  The  effective  chord  and 
centerline  of  the  hydrofoil  in  the  absence  of  injection  (but  in- 
cluding boundary- layer  effects!  are  shown  in  detail  in  Figure  l4. 
As  a first  approximation  it  is  appropriate  to  assume  that  the 
injection  alters  the  boundary-layer  thickness  only  on  the  side 
on  which  the  injection  is  made  and  only  downstream  from  the  point 
of  injection  . The  effective  centerline  of  the  hydrofoil  will 
then  be  displaced  vertically  by  an  amount  equal  to  A5  where 
A6V  is  the  change  (assumed  positive  when  there  is  thickening  on 
the  upper  surface  or  a thinning  on  the  lower  surface)  in  the 

* ’/or  the  sake  of  convenience  only  a general  thinning  Is  shewn  in 
Mgs.  I'l  and  1 *>.  As  mentioned  earlier,  there  can  be  a thinning 
or  thickening  of  the  boundary  layer  locally  depending  on  the 
relative  magnitudes  of  the  effects  due  to  the  adverse  pressure 
gradient  and  the  polymer. 
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boundary-layer-displacement  thickness  due  to  injection.  The 
effective  centerline  and  chord  are  shown  as  solid  lines  in 
Figure  14,  with  T"  being  the  new  effective  trailing  edge  of  the 
hydrofoil . 

It  can  be  seen  from  Figure  14  that  the  injection  leads  to 


a change  in  the  effective  angle  of  attack  by  an  amount  equal  to 
-A6*/2c,  where  A6*  is  the  change  in  boundary-layer- displacement 
thickness,  due  to  injection,  at  the  trailing  edge.  It  can  also 


be  readily  shown  that  the  change  in  camber  is  given  by  -l(A5*  - 
xA5£) . Since  the  camber  effect  depends  on  the  details  of  the 
interaction  between  the  boundary  layer  and  the  injected  flow,  it 
may  be  either  positive  or  negative  locally  and  change  sign  with 
downstream  distance  depending  on  the  specific  nature  of  the 
interaction.  As  already  pointed  out,  the  fact  that  different 
polymers  lead  to  significantly  different  results  implies  that 
the  characteristic  interaction  length  is  likely  to  be  dependent 
on  the  parameter  Vmt , where  t is  the  polymer  relaxation  time, 

V 

It  may  be  executed  that  the  nondimeneional  ratio  6..  /V**,  where 
6j.  is  the  displacement  thickness  at  the  location  of  the  injec- 
tion, will  be  an  important  parameter. 


Of  course,  the  exact  nature  of  the  interaction  of  the  injected 


polymer  solutions  with  the  boundary  layer  or  with  the  outside  flow 


Is  unknown,  and  it  is  not  the  objective  of  the  present  study  to 
attempt  to  describe  the  details  of  this  interaction.  It  is  our  oh  jot; 
tlv©  to  merely  Illustrate  that  such  Intiraotienk^ait'bc  the 
plausible  reasons  for  the  observed  changes in  ■thtf^pressure  ;dis-. 
tributions,  lift  and  drag,  Kvon  if  the  injected  flow  enters  t'ue 
outside  flow  as  a “swollen  jet”,  the  computational  procedure  de- 
scribed heroin  can  be  utilised  provided  that  the  displacement 
effect  of  this  jet  is  known. 
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It  should  be  noted  that  the  circulation  around  the  modi- 
fied centerline  LTn  will  not  be  that  specified  by  the  Kutta- 
Joukowski  hypothesis^  nor  will  it  be  merely  the  circulation  in 
the  absence  of  injection  (that  is,  for  the  centerline  LT')  cor- 
rected for  the  effects  of  the  changed  angle  of  attack  and  camber. 
The  new  value  of  the  circulation  can,  and  should  be,  uniquely 
determined  by  using  the  condition  that  the  velocities  at  the 
edge  of  the  (changed)  top  and  bottom  boundary  layers  at  the 
trailing  edge  must  be  the  same. 

Finally,  the  effect  of  the  injection  on  the  thickness  dis- 
tribution can  be  obtained  by  subtracting  A6*/2  from  the  basic 
distribution  when  there  is  a thinning  of  the  boundary  layer  and 
by  adding  it  when  there  .is  a thickening. 

The  four  separate  effects  can  be  summarized  as  follows : 

(i)  A change  in  angle  of  attack  due  to  a change  in 
the  boundary-layer  thickness  at  the  trailing 
edge)  given  by  Aa  = -A6{*  /2c. 

(ii)  A change  in  the  effective  camber  of  the  hydro- 
foil due  to  injection. 

(iii)  A change  in  circulation  duo  to  injection,  deter- 
mined by  using  the  theorem  stated  by  G.  I*  Taylor 
that  equal  and  opposite  amounts  of  vortlelty  must 

pass  into  the  wake  from  the  top  and  bottom  boundary 

, 11,17 

layers  * 1 . 

(iv)  A change  in  thickness  distribution  duo  to  injection. 

In  passing,  it  may  be  noted  that  in  their  analysis  of 

7 

boundary-layer  effects,  Fruman,  Tulin  and  Liu'  only  considered 
the  first  of  the  above  four  effects,  while  in  the  comparison 
shown  in  Figure  12  (that  is,  in  K qua t ion  [ 9 1 ) > all  of  the  ob- 
served lift  change  was  attributed  to  the  third  effect. 
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Before  making  quantitative  estimates  of  the  four  separate 
effects  of  injection,  it  is  first  relevant  to  consider  their 
qualitative  behavior  for  a specific  example,  say,  that  of  poly- 
mer injection  on  the  upper  surface.  If  the  effect  of  polymer 
injection  is  to  make  the  boundary  layer  downstream  of  the  point 
of  injection  behave  in  a more  laminar  fashion,  then  accompanying 
the  drag  reduction,  there  is  also  a reduction  in  the  boundary - 
layer-displacement  thickness  at  the  trailing  edge.  Thus,  there 
will  be  an  increase  in  the  effective  angle  of  attack.  Of  course, 
the  converse  is  true  if  polymer  injection  is  made  on  the  lower 
surface. 

As  already  pointed  out,  the  si gn  of  the  camber  effect  will 
depend  on  the  detailed  nature  of  the  Interactions  between  the 
Injected  and  boundary- layer  flows.  :in  general,  if  the  thinning 
effect  Increases  mono  tonic  ally  with  increasing,  downstream  dis- 
tance, as.  depicted  In  figure  l‘>(u),  then  there  will  be  a net 
positive  Increase  in  camber  (due  to  a decrease  In  the  negative 
camber)  as  shown.  However,  if  the  maximum  degree  of  thinning 
occurs,  at  a short  distance  downstream  of  the  point  of  Injection, 
as  shown  in  figure  l‘.>(b),  the  distribution  of  the  change  In  cam- 
ber may  be  negative  over  a large  portion  of  tie*  hydrofoil,  so 
that  the  net  effeet  Is  a reduction  in  lift.  If  tie?  maximum  de- 
gree of  thinning  oeeurs  at  a location  close  to  the  trailing  edge 
but  ahead  of  It,  as  shown  in  figure  th(e),  then  the  maximum 
(negative)  change  in  camber  will  also  occur  close  to  the  trail- 
ing edge,  gineo,  as  Is  well  known  (see,  for  example,  Ueferenee 
itf,  page  M),  the  values  of  the  camber  near  the  trailing  edge 
have  a more  profound  influer.ee  on  lift  than  the  valuer,  elsewhere, 
the  reduction  of  lift  in  this  case  is  likely  to  be  substantial . 

The  change  in  circulation  due  to  polymer  Injection,  say,  at 
the  top  surface,  can  be  determined  using  the  condition  that  the 
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new  value  of  the  circulation  should  be  such  that  the  velocity 
at  the  new  edge  of  the  top  boundary  layer  at  the  trailing  edge 
should  be  the  same  as  that  at  the  (unchanged)  edge  of  the  bottom 
boundary  layer  at  the  trailing  edge.  For  a hydrofoil  with  a' 
finite  trailing-edge  angle,  the  trailing  edge  is  a stagnation 
point,  so  that  near  the  trailing  edge  the  velocity  will  increase 
with  increasing  distances  from  it.  Thus,  if  the  circulation 
remained  unchanged  after  injection,  the  velocity  at  the  edge 
of  the  (thinned) top  boundary  layer  will  be  smaller  than  that 
at  the  bottom]  to  increase  this  velocity  to  equal  that  at  the 
bottom,  the  circulation  will  have  to  increase.  In  other  words, 
a thinning  of  the  top  boundary  layer  must  necessarily  be  accom- 
panied by  an  increase  in  the  effective  circulation  around  the 
hydrofoil* . 


Conversely,  a thinning  of  the  lower  boundary  layer  must 
necessarily  bo  accompanied  by  a decrease  in  circulation,  since 
to  make  the  velocity  at  the  new  edge  of  the  (thinner)  lower 
boundary  layer  at  the  trailing  edge  equal  to  that  at  the  (un- 
changed) upper  edge,  the  circulation  must  be  decreased. 


In  summary,  if  the  effect  of  polymer  injection  on  the 
boundary  layer  is  to  thin  it  downstream  of  the  point  of  injec- 
tion, the  angle-of-attack  effect  will  increase  the  lift  for 
upper  side  injection  and  reduce  it  for  lower  side  injection. 
The  camber  effect  may  either  increase  or  decrease  the  lift  de- 
pending on  the  details  of  the  interaction  between  the  injected 


■* Since  the  velocity  distribution  in 
edge  is  dependent  on  the  hydrofoil 
trailing  edge  angle  l.n  'particular, 
fluence  the  change. in  circulation', 
in  .Spence 1 s theory'0. 


the  region  of  the  trail  i.ng 
geometry  In  general,  and  the 
these  parameters  will  :l,n- 
Thoso  features  are  included 
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flow  and  the  boundary  layer.  The  change  in  lift  due  to  the 
modified  trail. ling- edge  condition  (Taylor's  theorem)  is  posi- 
tive for  upper-side  injection  and  negative  for  lower-side  in- 
jection. Thus  the  total  change  in  lift  may  be  either  positive 
or  negative  depending  on  the  relative  magnitude  of  the  differ- 
ent effects,  though,  in  general,  upper-side  injection  can  be 
expected  to  produce  an  increase  in  lift  and  lower-side  injec- 
tion a reduction  in  lift. 

It  is  now  appropriate  to  consider  some  quantitative  esti- 
mates. Equations  [3]  and  [4]  represent  the  induced  velocity 
distribution  for  a symmetric  hydrofoil  with  an  arbitrary  cir- 
culation r.  Eor  a cambered  hydrofoil  the  Induced  velocity  due 
to  thickness  distribution  is  still  represented  by  Equation  [4] 
However,  utilising  the  method  developed  by  Weber  * ^ it  can 
be  shown*  that  the  induced  velocity  duo  to  the  vortex  distrilm 
tion'for  a cambered  hydrofoil  is  given  by 


Where,  as  before,  the  positive  sign  applies  to  the  upper 
surface1  and  the  negative  sign  to  the  lower  surface.  Also, 
a |.,j.  :ls  the  angle  of  attack  measured  with  respect  to  the 
effective  chord  line,  1..T'  , w is  the  height,  of  the  effective 
camber  line  above  the  x axis,  while  s,  Is  the  thickness  dis- 

* ih 

tr.UmUen  of  the  hydrofoil  (including  the  boundary-layer  dis- 
placement thickness) , 


* Hoe  Appendix  A, 
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Equations  [l4]  and  [6]  together  describe  the  velocity  dis- 
tribution on  the  hydrofoil.  From  these  equations,  the  change, 
AV  , in  the  velocity  distribution  due  to  injection  can  be  ob- 
tained as. 


AV 
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where,  as  before,  the  positive  sign  applies  to  the  upper  sur- 
face and  the  negative  sign  to  the  lower  surface.  The  four 
terms  on  the  right-hand  side  represent  the  effects  on  the 
velocity  distribution  of  changes,  respectively,  in  the  effec- 
tive angle  of  attack,  in  circulation,  in  camber  and  in  the 
thickness  distribution.  Note  that  the  first  three  terms  are 
of  opposite  signs  on  the  two  sides  of  the  hydrofoil , while 
the  last  term  is  of  the  same  sign  on  both  sides.;  thus  the  last 
term  will  net  influence  the  lift  change. 


Even  though  Equation  [15]  appears  to  be  quite  complex, 
the  changes  in  velocity  can  be  calculated  fairly  readily  if 
the  changes  in  boundary -layer  displacement  thickness  due  to 
injections  are  known,  since  simple  algebraic  methods  for  evalu- 
ating the  Intergrals  appearing;  in  :lt  are  available  in  the  liter 
Hturo'L',,^>,'l'‘',  Equation  [15]  can  be  simplified  by  utilising 
the  expressions  given  earlier  relating  the  change  in  boundary 
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layer  thickness  to  changes  in  the  effective  angle  of  attack, 
camber  and  thickness  distribution.  After  some  algebra,  one  ob- 
tains. 


W=±I^  *L  + 


'/  VSi 


'l  + (dz  / dx): 


dA6* 


He1  1 dx 1 
1-x1  f x-x! 


[16] 


where  V and  V,  represent,  respectively,  the  velocities  on 
the  top  and  bottom  surface  of  the  hydrofoil.  In  the  deriva- 
tion of  this  equation  the  approximation  hwl  has  been  made  in 
the  second  term,  so  as  to  facilitate  physical  understanding 
of  the  phenomena  involved.;  the  approximation  will  not  signifi- 
cantly alter  the  numerical  values  of  the  velocity. 

The  first  term  on  the  right-hand  side  of  Equation  [l6] 
represents  an  antisymmetric  effect  of  boundary-layer  dis- 
placement, while  the  second  term  reflects  both  symmetric  and 
antisymmetric  effects.  The  first  term  accounts  for  the  change 
in  circulation  due  to  a changed  traillng-edge  condition,  and 
produces  velocity  changes  which  are  equal  in  magnitude  but 
opposite  in  sign  on  the  two  sides  of  the  hydrofoil.  The  lift 
effect  represented  in  this  term  is  essentially  the  same  as 
that  in  Equation  [9]  and,  as  already  noted,  ascribing  all  of 
the  observed  lift  changes  to  this  effect  results  in  a pre- 
dicted pressure  distribution  which  is  in  good  agreement  with 
the  measured  value,  except  for  the  presence  in  the  latter  of 
certain  peaks  located  near  the  point  of  maximum  hydrofoil  thick* 
ness  and  on  the  same  side  as  the  one  on  which  the  injection  is 
made . 
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While  the  first  term  on  the  right-hand  side  of  Equation 
[16]  is  not  of  a form  which  can  yield  the  experimentally  observed 
pressure  peaks,  the  second  term  is  indeed  of  this  form,  as  can 
be  seen  from  the  following  argument.  If  we  assume  that,  to  the 
first  order,  injection  changes  the  boundary-layer  thickness  only 
on  the  side  on  which  the  injection  is  made,  A6*  and  A6 * will  be 
equal  in  magnitude  as  well  as  in  sign  when  injection  is  on  the 
top  surface,  while  they  will  be  equal  in  magnitude  but  of  opposite 
. sign  when  injection  is  on  the  bottom  surface;  see  Equation  [12]. 
First  consider  tne  case  of  injection  of  polymer  or  water  on  to 
the  top  surface.  For  this  case  it  can  be  seen  from  Equation 
[l6]  that  the  value  of  the  integrand  in  the  second  term  has  a 
singularity  at  x = x',  while  it  is  bounded  at  all  values  of  x on 
the  lower  surface*.  Thus,  depending  on  the  specific  behavior  of 
A6*,  it  is  likely  that  a pressure  peak  of  the  type  observed  in 
the  experiments  can  occur  on  the  top  surface,  but  not  on  the 
bottom  surface.  Also,  assuming  that  polymer  and  water  injections., 
produce  changes  of  opposite  sign  in  the  displacement  thickness 
(an  assumption  which  is  consistent  with  the  observed  drag  behavior), 
it  can  be  seen  that  water  injection  will  produce  peaks  of  sign 
opposite  to  those  produced  by  polymer  injection. 

‘X-  -X- 

For  bottom  side  injection,  A6  and  A6  are  of  opposite  sign 
(though  of  equal  magnitude)  so  that  the  Integrand  in  Equation 
[l6]  is  bounded  everywhere  on  the  top  surface  but  not  on  the 
bottom  surface.  Hence  in  this  case  the  pressure  peaks  can  bo 
expected  to  occur  on  the  lower  surface  and  not  on  the  upper  sur- 
face, The  features  noted  above  are  exactly  the  ones  displayed 
by  the  experimental  data. 


*In  this  and  all  other  singular  integrals  occurring  in  this  ropor 
it  is  understood  that  the  Cauchy  principal  value  is  the  one  of 


t 


> 


concern. 
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Since  the  detailed  nature  of  the  interaction  of  the  injected 
flow  with  the  boundary  layer  is  not  known,  Equations  [153  and 
Cl6]  cannot  be  used  to  make  quantitative  predictions  of  the  veloc-- 
ity  changes  that  will  occur  due  to  injection.  However,  these 
equations  can  be  used,  along  with  the  measured  pressure  changes, 
to  estimate  the  shape  of  the  effective  hydrofoil  which  corresponds 
to  the  observed  changes-,  as  follows: 

The  change  in  the  effective  angle  of  attack  is  first  calcu- 
lated using  Equation  [13]  and  the  momentum  relation  expressing 

the  change  in  boundary  layer  thickness  at  the  trailing  edge  in 

7 

terms  of  the  change  in  skin-friction  drag  . This  last  quantity 
is  itself  determined  by  computing  the  change  in  the  pressure  drag 
from  the  measured  pressure  distributions  and  by  subtracting  this 
from  the  measured  changes  in  the  total  drag.  The  first  term  on 
the  right-hand  side  of  Equation  [ 153  can  then  be  readily  calculated. 


If  the  total  velocity  field  around  the  hydrofoil,  including 
the  .boundary-layer-displacement  effect,  were  known,  then  the  change 
in  circulation  Ah  can  be  calculated  using  the  trailing-edge  condi- 
tion that  the  velocities  at  the  edges  of  the  top  and  bottom  bound- 
ary layer  be  the  same*.  Lacking  such  information,  the  second  term 
on  the  right-hand  side  of  Equation  [15],  which  is  the  same  as  the 
first  term  in  Equation  [16],  can  be  estimated  for  various  assumed 
values  of  Ah.  The  fourth  term  on  the  right-hand  side  of  Equation 
[153  can  be  eliminated  by  taking  the  differene e between  the  mea- 
sured velocity  changes  on  the  top  and  bottom  surfaces  of  the. hydro- 
foil. Thus,  finally,  the  third  term  (which  represents  the  camber 
change)  on  the  right-hand  side  of  Equation  [1 6]  can  be  expressed 


in  terms  of  the  measured  pressure  changes.  The  method  of  Weber 


.15,16 


*The  method  of  Spcncol3  provides  an  approximate  procedure  for  carry- 
ing out  the  step.  However,  oven  this  approximate  procedure  la  too 
complex  to  be  useful  here, 
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can  then  he  used  to  calculate  the  camber,  since  she  has  tabulated 
the  influence  coefficients  necessary  for  the  inversions  of  the 
integral. 

Figure  l6  shows  the  results  of  three  sample  calculations  for 
the  cases  of  Jaguar,  Polyacrilamide  and  Polyox  injections  on  the 
suction  side  of  the  hydrofoil  at  10 # chord,  2.5°  angle  of  attack 
and  at  an  injection  velocity  equal  to  thirty  percent  of  free-stream 
velocity.  Calculations  are  shown  both  for  Ah  - 0 and  All  = 0,015. 

As  expected,  the  calculations  show  that  significant  changes  in 
camber  occur  only  downstream  of  the  point  of  injection.  For  both 
of  the  assumed  values  of  Ah  in  Figure  l6.  Jaguar  injection  pro- 
duces the  maximum  positive  camber,  while  Polyacrilamide  and  Polyox 
injections  give  successively  lower  values.  It  is  relevant  to  note 
that  the  change  in  camber  is  not  confined  to  a small  region  close 
to  the  point  of  injection,  but. occurs  over  the  entire  downstream 
region  of  the  hydrofoil.  The  changes  with  downstream  distance  in 
sign  of  the  camber  (especially  evident  in  the  cases  for  Ah  - 0.015) 
may  be  due  to  the  competing  influences  of  the  adverse  pressure 
gradient  and  polymer  relaxation. 

Also  shown  in  Figure  16,  for  the  sake  of  comparison,  are  the 
results  for  the  cases  of  30#  chord  suction- side  and  10#  chord 
pressure-side  Injections  of  Jaguar  at  V^/VM  --  .3  rate  and 
a - 3.25°.  As  expected,  the  changes  in  camber  for  the  case  of 
pressure-side  injection  are  negative  (corresponding  to  the  observed 
decrease  in  lift).  Xt  is  interesting  to  note  that  for  the  case  of 
thirty-percent-chord  injection,  substantial  changes  in  camber  begin 
somewhat  ahead  of  the  point  of  injection. 

The  sample  calculations  discussed  above  illustrate  that  the 
changes,  due  to  injection,  in  the  cumber  corresponding  to  the 
observed  changes  in  the  pressure  distribution,  lift  and  drug  are 
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quite  complex.  Thus,  even  when  the  effects  of  polymer  injection 
are  viewed  in  terms  of  the  changes  in  the  "effective"  shape  of  the 
hydrofoil,  no  straightforward  phenomenological^  explanation  suggests 
itself.  However, 'the  observed  changes  in  pressure  distribution  and 
the  corresponding  calculated  changes  in  camber  are  consistent  with 
those  that  can  be  expected  from  a complex  interaction  of  the  in- 
jected flow  with  the  boundary'  layer  and  the  external  pressure 
gradient.  : \ ' v ' 


V ■ t 
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IV. 2 Correlations  of  Lift  and  Drag  Data 

As  already  pointed  out,  the  fact  that  under  otherwise  identi- 
cal conditions  different  polymers  lead  to  different  lift  and  drag 
effects  implies  that  viscoelastic  effects  may  be  responsible,  at 
•least  in  part,  for  the  observed  behavior.  Based  on  similar  corre- 
lations that  exist  in  Pitot  tube  data  for  dilute  polymer  flows, 
Fruman  and  Tulin1^  have  suggested  that  the  lift  and  drag  effects 
should  be  correlatable  in  terms  of  a nondimensional  parameter 
formed  by  dividing  the  relaxation  length  V^t  by  a characteristic 
geometric  length  scale  in  the  problem,  where  V±  is  the  polymer  in- 
jection velocity  and  t is  the  relaxation  time.  Such  a correlation 
is  shown  in  Figure  17..  where  all  of  the  lift  data  are  shown  plotted 
against  the  parameter  (V^t/c),  with  c being  the  hydrofoil  chord. 

Note  that  the  ordinates  in  Figure  17  represent  the  difference 
in  the  lift  changes  between  polymer  and  water  injections j this  is 
a logical  way  of  correlating  the  data  since  the  relaxation  time 
corresponding  to  water  injection  is  zero  and  since  for  small  values 
of  the  abscissa,  the  ordinate  can  be  expected  to  tend  to  aero. 


While  there  is  considerable  data  scatter  in  Figure  17,  it  is 
nevertheless  clear  that  the  correlation  :ls  a valid  one  since  all 
of  the  different  curves  shown  have  essentially  the  same  shape.  It 
is  relevant  to  note  that  each  of  the  curves  given  in  Figure  17  not 
only  correlates  the  cases  for  the  different  polymer  injections, 
but  also  the  different  injection  velocities.  The  corrcla tilons 
also  show  that,  in  general,  at  a fixed  value  of  the  abscissa  larger 
hydrofoil  angles  display  larger  lift  changes,  and  that  the  effects 
are  somewhat  smaller  for  the  30*;'*  chord  injection  than  for  the  10# 
chord  injection. 


The  data  show  that  for  either  injection  position  the  maximum 
lift  effects  occur  when  the  value  of  the  abscissa  is  between  0.03 
and  0.04,  with  suction-side  injection  giving  a lift  increase  and 
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pressure-side  injection  a lift  decrease.  As  the  value  of  the 
abscissa  increases  beyond  this  critical  value,  the  absolute  magni 
tude  of  the  lift  effect  decreases,  and  the  sign  of  the  effect 
actually  reverses  'beyond  a certain  value  of  the  abscissa.  Thus, 
the  correlation  seems  to  offer  an  empirical  explanation  as  to  why 
and  how  the  exception  to  the  rule  that  "suction-side  injection 
increases  lift  and  pressure-side  injection  reduces  lift"  arises. 


It  was  mentioned  earlier  that  Fruman,  Tulin  and  Liu^  were 
able  to  successfully  correlate  the  observed  changes  in  lift  for 
different  free-stream  velocities  by  plotting  the  observed  lift 
force  against  the  local  velocity  at  the  injection  slit.  Thus  it 
is  logical  to  attempt  to  correlate  the  data  in  the  present  case 
by  normalizing  the  observed  lift  changes  not  by  the  free-stream 
velocity,  but  by  the  local  velocity  at  the  location  of  the  injec- 
tion slit.  Such  a correlation  is  shown  in  Figure  l8.  Again, 
while  there  is  considerable  data  scatter,  the  trends  in  the  data 
are  readily  apparent.  It  is  noteworthy  that  all  of  the  four 
curves  shown  have  similar  shapes.  The  feature  that,  in  general, 
the  30^-chord  injection  yields  smaller  lift  effects  than  the  10;j- 
chord  Injection,  may  be  related  to  the  fact  that  in  the  former 
case  less  of  the  hydrofoil  surface  is  affected  by  the  injection. 


However,  a further  correlation  of  the  data  to  account  for  the 
the  injection  position  is  not  readily  apparent. 


A correlation  for  the  relative  drag  change  due  to  Injection, 
similar  to  that  shown  in  Figure  17  for  lift,  is  shown  in  Figure 
19  in  terms  of  the  parameter  V^V'e.  The  shape  of  the  curves  are 
seen  to  be  the  same  as  those  for  the  lift  effect,  with  the  occur- 
rence of  the  maximum  drag  reduction  coinciding  with  the  maximum 
lift  reduction  for  pressure-side  Injection,  and  the  occurrence  of 
the  minimum  drag  reduction  (or  maximum  drag  increase)  coinciding 
with  the  maximum  lift  increase  for  suction-side  injection. 


, . \ 
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Tii  e changes  in  drag  due  to  injection  can  be  resolved  into 
two  components,  namely,  those  due  to  changes  in  the  pressure  and 
frictional  parts  of  the  total  drag  of  the  hydrofoil.  The  change 
in  the  pressure  drag  -can  be  computed  from  the  measured  changes  in 
the  pressure  distribution  and  from  the  known  shape  of  the  hydrofoil 
by  taking  the  component  of  the  local  pressure  force  in  the  direction 
of  the  free  stream.  The  change  in  the  friction  drag  can  then  be 
deduced  by  subtracting  the  calculated  pressure  drag  change  from 
the  measured  change  in  the  total  drag. 


A resolution  of  the  drag  components  is  shown  in  l-’igure  JO 
for  the  case  of  Polyox  injection  at  ten  percent  of  free- stream 
velocity.  It  can  be  seen  that  for  the  lower  injection  velocity 
the  pressure  drag  actually  increases  both  for  suction-side  and 
pressure-side  injections  at  all  values  of  the  angle  of  attack. 

The  decrease  in  skin-friction  drag  seems  to  bo  nearly  the  same 
for  inj oo lions  on  either  side  of  the  hydrofoil  and  for  all  angles 
of  attack.  It  is  interesting  to  note  that  the  fact  that  the 
total  drag  reduction  for  suction-side  injection  at  is  small 
(see  Figure  3)  is  because  the  increase  in  pressure  drag  nearly 
offsets  the  reduction  In  skin-friction  drag. 

For  injection  at  the  3df'-chord  location,  the  change  in  pres- 
sure drag  is  negative  at  small  angles  of  attack  for  -injection  on 
either  side  of  the  foil,  though  It  becomes  positive  at  larger 
angles  of  attack.  The  reduction  in  skin-friction  drag  In  this 
ease  is  somewhat  smaller  than  that  in  the  10#- chord  injection 
ease,  consistent  with  the  fact  that  a smaller  fraction  of  the 
hydrofoil  surface. area  is  affected  in  this  ease. 

Figures  Id  and  .'0  also  Illustrate  that  the  drag  and  lift 
effects  are  closely  related.  Indeed,  If  the  pressure  drag  com- 
ponent alone  Is  correlated  with  (\h  r/o) , It  is  clear  that  cor- 
relations of  the  type  shown  in  Figures  17  and  K>  must  result. 
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V.  CONCLUDING  REMARKS 


As  is  evident  from  the  discussion  given  in  the  previous 
sections,  the  effects  of  polymer  injections  on  lifting  surfaces 
involve  fairly  complex  phenomena.  While  the  actual  magnitude 
of  the  lift  changes  are  in  themselves  not  large,  it  is  neverthe- 
less necessary  to  seek  a rational  explanation  for  the  observed 
behavior  of  the  lift-drag  ratio  before  the  use  of  polymer  addi- 
tives to  enhance  the  performance  of  propellers  and  hydrofoil 
craft  can  become  a reality. 


In  the  present  study,  the  observed  lift  effects  were  viewed 
as  arising  from  two  separate  steps  of  interaction  of  the  injec- 
ted flow  with  the  flow  about  the  hydrofoil.  In  the  first  stop 
the  injected  flow  was  assumed  to  interact  with  the  boundary  lays 
and  the  external  pressure  gradient  in  some  (as  yet)  unspecified 


manner  and  change  the  "effective"  shape  of  the  hydrofoil.  This 
step  of  the  interaction  is  undoubtedly  dependent  on  the  visco- 
elastic characteristics  of  the  injected  polymer  since  under 
otherwise  identical  conditions,  different  polymers  give  consider 
ably  different  results;  a detailed  consideration  of  this  inter- 
action was  outside  the  scope  of  the  present  study. 

In  the  second  step  of  the  interaction,  the  changed  " effect! 
hydrofoil  shape  was  assumed  to  produce  changes  in  the  pressure  a: 
velocity  distributions  on  the  hydrofoil.  It  was  also  assumed  th* 
the  latter  changes  can  be  calculated  utilising  potential  theory 
without  any  further  considerations  of  the  viscoelastic  effects. 
In  terms  of  the  ability  to  theoretically  predict  the  exact  nutur* 
of  the  lift  and  pressure  changes  that  can  occur  due  to  injection. 


solution  of  the  first  step  is  a necessary  prerequisite  to  the 
utilisation  of  the  second  step,  and  this  was  outside  the  scope  o 
the  present  study.  However,  the  measured  pressure  distributions 
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as  well  as  the  observed  lift  and  drag  changes  can  be  used  along 
with  potential  theory  to  calculate  the  "effective"  hydrofoil 
shape  that  would  produce  the  type  of  changes  measured.  In  turn, 
the  effective  shape  can  be  examined  to  see  if  it  is  compatible 
with  what  can  be  expected  based  on  plausible  hypotheses  of  inter- 
actions between  the  injected  flow  and  the  boundary  layer  as  well 
as  the  external  pressure  gradient,  without  having  to  consider 
the  details  of  the  interactions  themselves.  This  is  the  approach 
that  was  utilized  in  the  present  study. 

Analysis  of  the  types  of  interactions  that  can  take  place 
between  the  injected  flow  and  the  external  flow  around  the  hydro- 
foil indicates  that  four  different  types  of  changes  in  the  pres- 
sure distribution  can  occur.  These  are:  a change  in  the  effec- 
tive angle  of  attack  due  to  a change  in  the  boundary-layer  dis- 
placement thickness  on  the  side  on  which  the  injection. 1 s made; 
a change  in  the  effective  camber  of  the  hydrofoil  due  to  a change 
in  the  boundary-layer  displacement  thickness;  a change  in  the 
circulation  due  to  a change  in  the  trailing-odge  condition  re- 
quiring that  the  velocities  at  the  edges  of  the  top  and  bottom 
boundary  layers  at  the  trailing  edge  must  be  the  same;  and,  a 
change  in  the  effective  thickness  distribution  of  the  hydrofoil . 

The  qualitative'  features  of  the  observed  pressure  distribu- 
tions as  well  as  the  lift,  and  drag;  changes  art'  compatible  with 
the  types  of  Interaction  described  above.  Calculations  carried 
out  using  potential  theory  to  determine  the  "effective”  shape  el* 
the  hydrofoil  that  produces  the  observed  changes  are  also  com- 
patible with  the  types  of  interaction  described.  However,  lack- 
ing; a specific  theory  for  the  interactions  themselves,  the  ex- 
planation offered  herein  has  to  bo  regarded  as  tentative  only. 

On  the  other  hand,  the  following;  points  support  the  present  ex- 
planation, and  these  have  to  bo  taken  into  account  in  the  develop- 
ment of  any  alternative  explanations. 
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(i)  Even  though  in  many  cases  polymer  injection  on  the 
suction  side  of  the  hydrofoil  produces  a lift  increase  and  that 
on  the  pressure  side  a lift  decrease,  this  is  not  always  the 
case.  In  general,  for  injection  on  either  side  of  the  hydro- 
foil the  lift  can  increase  or  decrease  depending  on  the  polymer, 
the  injection  velocity,  the  location  of  the  injection,  and  the 
angle  of  attack. 

(ii)  Both  polymer  and  water  injections  produce  peaks  (of 
opposite  signs)  in  the  pressure  distribution  on  the  side  on 
which  the  injection  is  made.  These  peaks  do  not  occur  immedi- 
ately downstream  of  the  injection  slit,  but  occur  near  the 
location  of  the  maximum  hydrofoil  thickness  for  all  cases. 

(iii)  Pressure  distributions  in  the  absence  of  injection 
indicate  the  presence  of  boundary-layer  effect  on  the  circula- 
tion around  the  hydrofoil. 
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FIGURE  3 - EFFECT  OF  ICf;  CHORD  INJECTION  OF  V/ATER  AND  POLYMERS 
ON  THE  DRAG  COEFFICIENT  OF  THE  HYDROFOIL 
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FIGURE  5 - EFFECT  OF  10#  CHORD  INJECTION  OF  WATER  AND  POLYMERS 
ON  THE  LIFT  COEFFICIENT  OF  THE  HYDROFOIL 


FIGURE  6-  EFFECT  OF  30 2 CHORD  INJECTION  OF  WATER  AND  POLYMERS  ON  THE 
LIFT  COEFFICIENT  OF  THE  HYDROFOIL 
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FIGURE  7 - DIFFERENCE  IN  PRESSURE  COEFFICIENT,  AC  VERSUS 
CHORD  DISTANCE,  x /c  p 
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FIGURE  9 - COMPARISON  BETWEEN  MEASURED  AND  CALCULATED  PRESSURE 
DISTRIBUTIONS 


FIGURE  10  - COMPARISON  BETWEEN  MEASURED  AND  CALCULATED 
PRESSURE  COEFFICIENTS  ON  THE  UPPER  AND  LOWER 
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(B)  WATER,  V./V^  .3 
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FIGURE  12-  COMPARISON  OF  MEASURED  CHANGES  IN  PRESSURE  Will! 
THOSE  CALCULATED  FOR  A SYMMETRIC  HYDROFOIL 
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(a  j SCHEMATIC  OF  EFFECT 


(b)  CAMBER  EFFECT 


FIGURE  13  - BOUNDARY-LAYER  DISPLACEMENT  EFFECT  ON  A 
SYMMETPIC  HYDROFOIL 
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FIGURE  14  - EFFECT  OF  POLYMER  IN  JECTION  ON  BOUNDARY 
LAYER  DISPLACEMENT 
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(o  ) THINNING  INCREASES  MONCTON ICALLY  WITH  DOWNSTREAM  DISTANCE 


( b ) MAXIMUM  THINNING  OCCURS  CLOSE  TO  POINT  OF  INJECTION 


( c ) MAXIMUM  THINNING  OCCURS  CLOSE  TO,  BUT  AHEAD  OF,  TRAILING  EDGE 


FIGURE  15  - TYPES  OF  OOUNDARY  LAYER  INTERACTION 
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FIGURE  16  - CALCULATED  CAMBER  DISTRIBUTIONS 
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FIGURE  18  - LIFT  CHANGE  VERSUS 
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APPENDIX  A 


VELOCITY  DISTRIBUTION  INTEGRALS  FOR  SYMMETRIC 
AND  CAMBERED  HYDROFOILS 


IB 

Weber  v has  deduced  an  integral  relation  for  the  velocity 
distribution  on  a flat  plate  at  an  angle  of  attach,  a , by  using 
conformal  transformation.  In  particular,  she  has  utilized  the 
transformation 


t 35  1 1 


R2 


to  transform  the  flat  plate  - s £ £ 4-  into  a circle  (for  the 
notation,  see  Reference  15).  The  velocity  distribution  on  one  f. 
plate  for  any  arbitrary  circulation  F around  it  is  then  exp res sib 
in  terms  of  the  complex  velocity  potential,  F , as,  . 


For  the  Ku  1 1 a - J oukow o hi  condition  to  be  satisfied, 
edge,  given  by  ; • i , must  be  a stagnation  point,  thus 
that  T = rlc  = IttV^  a.  Then,  in  the  physical  coordinate 
Equation  [A--]  gives, 

V = ' V.  • 


*•  *•’*  A *’  W>  * f i'  » 
^ • w , . u ^ ■ '-*«g 

X »>  *• 

» k^V*,  ■ * » k'  ,t 

***•■?  ^ M *« 

>3  evi.ti 


Uh 


For  an  arbitrary  circulation  P defined  by  Equations 
Equation  Ca-23  yields, 


VaiV„a 


h - x 

VmTT“x) 


r . 


1 


Of  course,  Equation  CA-4]  reduces  to  Equation  [A- 3.1  when  U « 1. 

For  a cambered  hydrofoil,  Weber1**  has  shown  that  the  velocity 
distribution  is  given  by, 

v = i|  • 
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when  the  circulation  is  specified  by  the  hutta- Joukowski  condition. 
When  the  circulation  is  specified  by  Equations  [1]  and  [2],  Equa- 
tion [A-5]  becor.es, 


V = + =•  V a 

TT  « 


h - x 
x(l  - x) 


1 dr  'x*  ' 
c ' 


- X1  X - x» 


. Ca-6] 


Note  that  Equations  [A-51  and  [A -6]  reduce,  respectively,  to 
.Equaticns  [A- 3]  and  [A-^2  when  is, /dx  = - a. 
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APPENDIX  B 


LIFT,  DRAG  AND  PRESSURE  MEASUREMENT  DATA 
(DECEMBER  8-19,  1975)  ON  THE  NACA 
634-020  TWO-DIMENSIONAL  HYDROFOIL 


VF  = free  stream  velocity  (Vffl  used  in  the  text) 
CD  = drag  coefficient 

CL  - lift  coefficient 

X/C  = chordwise  pressure-tap  position 
CPP  a pressure  coefficient  on  the  pressure  s3de 

CPS  = pressure  coefficient  on  the  suction  side 

Ending  N means  "no- injection"  case 
Ending  W means  "with-injeetion"  case 
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